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Abstract. We obtain cubature formulas of volume potentials over bounded domains com- 
bining the basis functions introduced in the theory of approximate approximations with their 
integration over the tangential-halfspace. Then the computation is reduced to the quadrature 
of one dimensional integrals over the halfline. We conclude the paper providing numerical tests 
which show that these formulas give very accurate approximations and confirm the predicted 
order of convergence. 

1 Introduction 

We consider the volume potential of modified Helmholtz operators in M. n 



with A G C. If A 2 ^ 0, then the fundamental solution K\(x) is given by 



where A £ C \ (- oo,0] and K v is the modified Bessel function of the second kind. The funda- 
mental solutions of the Laplacian (A = 0) are well-known 



A n = -A + A 2 




For /eC^fiC R n , the volume potential 

u(x) = /C A /(x) = / k a (x - y)/(y)dy (1.1) 
Jn 

provides a solution of 

f /(x) x g n 

otherwise 

We study cubature formulas for the volume potential (1.1) using the concept approximate ap- 
proximations (see [4]). The cubature of volume potentials over the full space and over high- 
dimensional halfspaces has already been studied in [1] and [2], respectively. In [5] cubature 
formulas based on approximate approximations for the single layer harmonic potential were 
considered. 

Assume / G C N (fl). We extend / outside Q with preserved smoothness and we denote by 
/ G C^(IR n ) the continuation of /. Assume that there exists C > such that 

ll/ilw^Rn) < C\\f\\ W N {n) . 

We introduce a uniform grid {/im} with step h. A cubature formula for (1.1) can be obtained 
if we replace / by the approximate quasi-interpolant 

Mh,vJ(x) = £ f(hm)r, (1-2) 

where rj G 5(M n ) and satisfies the moment conditions of order N 



/ 

JR 7 



r/(x) x Q cbc = 5 , a , < \a\ < N. (1.3) 
The quasi-interpolant (1.2) approximates / in U. It is known ([4]) that 

JV-l 

\f(x)-M h ,vf(*)\<c(Vvh) N \\X7 N f\\ LO c + Y,Zk(Vvh) k \X7 k f(x)\ 

k=0 

with 

s k < V \V k Tr)(Vv m )\; lim V |V fc Jr,(^m) = 0. 

meZ™\{0} m£Z"\{0} 

Since 77 is a smooth and rapidly decaying function, for any error e > one can fix r > and the 
parameter T> > such that the quasi-interpolant with nodes in a neighborhood of £1 



(1.4) 



d(hm,fl)<r hVv 
approximates / with 

|/(x) - Ml v fW\ = 0({Vvh) N + e)\\f\\ w s (1.5) 



2 



for all 

We use the quasi- interpolant (1.4) of the density to obtain a cubature formula of the volume 
potential (1.1) 

/C M /(x) = K x (Ml v f){x) = /(^ m ) j n ^ ~ y) 1 (^fif) ^ (L6) 

d(hm,£T)<r hy/V 

Since K\ is a bounded mapping between suitable function spaces, the differences K.\^f{x) — 
/Ca/(x) behave like estimate (1.5). We use the radial generating functions 

r ?2 M(x) = ^ 2 LK)(|x| 2 )e-W 2 , (1.7) 
where are the generalized Laguerre polynomials 

which satisfy the moment condition (1.3) with iV = 2M. Formula (1.6) would give a cubature 
of (1.1), if the integrals 

I KA(X - y)?? (^) dy 

could be computed efficiently for nodes with d(hm, Q) < r h\/V. 

For smooth domains, we propose to replace the integrals in (1.6) by integrals over the 
tangential-halfspace at a point of d£l with minimal distance to hm. It is proved that these 
formulas approximate (1.1) with the order O {(h\/V) 2 ) . We conclude the paper providing nu- 
merical tests which show that these formulas give very accurate approximations and the order 
of convergence cannot be improved. 



2 Cubature based on (1.2) 

In this section we study the approximation of the integral (1.1) over a bounded region $7 C W 1 
with smooth boundary by the sum (1.6) for appropriately chosen r > 0. 

Denote by Phm a point of d£l with minimal distance to hm and by np hm the normal at Phm 
directed towards the interior of Q. Let Th m be the halfspace bounded by the tangential plane 
at the point Phm such that the inner normal at Phm coincides with np hm . 

We define the following cubature formula for the volume potential (1.1) 

K M /(x) = 2^/ 2 £ f(hm) [ K A (x-y)^fc^W (2.1) 

To derive explicit expressions for the integrals in (2.1), as basis functions we use (1.7), satisfying 
the moment condition (1.3) with N = 2M (cf. [4]). We have 
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Lemma 2.1. Let Re A 2 > 0, n > 1 or Re A 2 = 0, n > 3. T/ie solution of the equation 

??2m(x), x n > a, 



(-A + X z )u= { 
is given by the one- dimensional integral 



o, 



x n < a 



F 2 (t,x n ,a) 



8vr™/ 2 J 
where 



e -A^/4 e -|x|V(i+t) (erfc(F(t,x TO ,a))P M (x,t) + e Q M (x,t,a)) rft (2.2) 



Af-l ^ 
P M (x,t) = ^2 (J + t )fc+n/2 



Q M (x,t,a) = 



2 

(i + t) (n-i)/2 2^^ {k _ iyA k- 



k 

M-l k 



((n-3)/2) 



2(k-l) 



((2{k-l)\ ( x n \ i?j_i(F(t,x n ,a)) /a- x n \H 2 ( k -l)-j(aT 



wi/i i/ie function 



F(t,x,a) = ] f^(a- T ^ t ). 
Hk denote the Hermite polynomials 



(2.3) 



(2.4) 



Proof. We consider the heat equation in M n 

<9 t z - Az = , z(x, 0) = < 
whose solution is given by the Poisson integral 



%m(x), x n > a, 
0, x n < a, 



z(x, i) 



(4yrf) n / 2 



-|x-y|V(4t) 



V2M(y) dy. 



Vn>a 



The representation ([4, p. 55]) 



M-l 



r ?2M (x) = vr-«/ 2 ^^A^e-N 2 



i=o 



(2.5) 



shows that 



M-l 



(-1) 



r72M(x',x n ) =7T n/2 ^ 

j /■-(> 



j\43 I )^' e ^2(j-l) 



dxn 
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Hence 

M-\ , j 



y > ( 4t )n/2^n z_> jW \ i r 3 iy ' y^.! y 



j=0 Z=0 



with 



From 



e ~ (x ~ y) ft -^2k e ~ y d y- ( 2 - 6 ) 

p 

f e -|x'-yT/(4t) A i e -|yr dy / = [ e -|x'-y'| 2 /(4t) e -|y'l 2 dy / 

= / _47r^x(n-l)/2 | x ,| 2/(1+4t) 

\l + 4tJ x 

and the relation 

A J e -\*'\ 2 /(W) - HOjilf e -|xT/(W) L ((n-3)/2) /JX? \ 

(i + t}» J vi + t/ 

([4, p. 121]), we obtain 



1 e -|x'| 2 /(l+4t) (_-Q 



(4t)n/2 ?r n (! + 4t )(n-l)/2 A> jl^ 
«=0 



From 



oo 

m (x,t,p) = J e-^^e-y 2 d y = ^l^-L-e-* 2 ^ erfc (F{t, x,p)) , 
p 

for /c > 1, integration by parts leads to 



,.2 



An i2i\ti £ 1 



y=p 



and the definition (2.4) gives 

d2W e -y 2 = (-ir-^e~y 2 H 2k t Ay) *- = ^ ^""^ H.( ^ 

In view of 

^_ erfc ( x ) = _L(_i)^ e - 2 ^ l(:r) ) ) £>i, 
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one gets for £ < 2k 



^— e erfc(F(t,x,p)) - 



d 



2k-£ 



(t(l + t)) k ~ e / 2 [dz 2k - e 
2 e - F2 ( t < x >p) 



erfc(z) 



^{t{l + t)) k -V 2 



z=F(t,x,p) 

H 2k -e-i{F(t,x,p)) . 



Therefore, since 



we obtain 

Q2k 



dx 2k 



ip (x,t,p) 



dx e (i + ty/ 2 



H 2 k 



2 (i + t)*+i/2 ' K \^T+t 



erfc(F(t, x,p)) 



^/i e -* 2 /(W) e -F\t,x, P ) 2 t^ (2k , (_!) 



(1 + t) k + 1 / 2 
Thus simple transformations give 



EC 



ip k (x,t,p) = e 



-x 2 /(l+t) 



erfc(F( t ,x, P ))g 2fc (-^=) 2(l + ; )fc+1/2 
^ (-1)* f(2k\ ( x x ^_!(F(t,x,p)) (P- x \tt ( \ 



(2.7) 



+ e -F 2 (t,x,p) 

Using (2.7), the relations 



H 2j (x) = {-\y^j\LY l,2 \x 2 ) and L f +b+1 \x + y) = ^L^L^y) 



(o+6+l) , 



1=0 



we find 



M-l 



E 



- j!4^ (l + ^+"/ 2 ^M 



E(|)(- 1 ) , «^r s,/2, (^)^-o( 7 fe) 



M-l 



- E 7TT^E^ , "- 3,/2, (^)4:! /2, (i £ fi) 



(1 + t)i+«/2 



leading to 



-|x| 2 /(l+4t) 

2tt™/ 2 



F 2 (4t,:r n ,a) 



erfc (F(4t, x n , a))P M (x, 4t) + -j= - Q M (x, 4t, a)) . 



□ 
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In the particular case n = 2 



I „ , N 1 / I Ixl 2 



P3 (x, t )=P 2 (x, f) + ^(^-2^ + J^) 

1 / 1 |x| 2 1 |x| 2 |x| 4 \ 

i Tt v 1 + r+t " (ITT) 2 + (TT*)2 " 2 (i + t)3 + 2(i + ty) ; 

Qi(x,t,a) = 0; Q 2 (x, t, a) = - - + ^ 3/2 (a + ^) 



1 Vt ( -2at ( x 2 \ /4|x| 2 -2a;| 7 2 \\ 

& < x ' t - "» = I TiTW { (TT7) + (» + TTl) (^17^ "TTi + 2a 5 )J- 

Remark 2.1. For sufficiently large \a\ > r the integrands in (2.2) are approximated by 



, a > r , 

-12/ 



2e -A 2 t/4 e -|x|V(i+t) p M ( X; t ) ) a < _ r , 
u>ii/i the error 0(e~ a ), which is in accordance with 

[ K X (x-y) V2M (y)dy = -^ f°° e'^ 4 e'^I^VM^ t) dt. 

Jr« 47T n / 2 Jo 

Theorem 2.1. Suppose that the generating function rj G 5(M n ) satisfies the moment condition 
(1.3) with N > 2. T/ien i/ie integral (1.1) is approximated by the sum 

jc XA f(x) = v-^ K hm )j KA(x - y) ^fe^) dy (2 - 8) 

with the error estimate 

|£ M /(x) - /C A /(x)| = 0(e + ci (/^) 2 ) (2.9) 

provided dfl has C 2 — smoothness. The saturation term e can be negligible small if V is large 
enough. 

Proof. We study the difference 

,„ < — , ~ P /v — hm \ 

dy 



1Cx,hfW ~ /CA, fc /(x) = P-"/ 2 £ /(H / « A (x - y) V i 1 -^) 

d ( hm ,n)<rhVv Jt ^ V hVT) J 



d(hm,Q)<rhVV 

as /l — >■ 0. 
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Since r?(y) is supported in the small neighborhood of the origin |y| < r, if hm G Q, such that 
d(hm, dft) > rh\[V we replace each integral 



/ K A (x-y)?7 
Jn 



y — hm 

hVV 



dy 



by integrals over r^ m without loss accuracy. 

Assume m G Z n : d(hm, d£l) < rh\[V. We choose a local coordinate system with the origin 
at the point Jim and the normal at the nearest point of 30 from coincides with the x n — 
axis (0', 1), 0' G R n_1 . In these coordinates system the halfspace is defined by T = T/j m = {y = 
(y',y n ) ■ y' € IR™" 1 ,?/™ > ph m } where 



P = Phm 



dist(0',dQ) ifO^n 
-dist(0',dn) ifOefi. 



We have \p\ < rh\fV. We assume that in a neighborhood U of the point (0', p) G dVL the domain 
0, is given by 

vn > f(y'), <p(o) = P , v^(o) = o. (2.io) 

We introduce the change of variable z = z(y) defined as 

z' = y', z n = y n - ip{y') + p 

and its inverse y = y(z) 

y' = z', y n = z n + ip(z') - p. 

Choose h such that B hr ^nQ C U and assume ro > r such that z(On5 ftrv ^) = Tf)B hrQ ^. 
Then the difference 

L ka((x - km) - y) * fe) dy - L ka((x - km) - y) * G^) dy 



takes the form 



/ 

JTnB 



ka((x — lim) — z) 77 



rnB .^(0) 



ka((x - /im) - z)) 



K A ((x-/im) -y(z))r/ 



f y(') \ 



y(g) 



dz 



+ 



/ [«a((x - /im) - z) - k a ((x - /an) - y(z)] 77 (7^= ) dz. 

J TnB rohV „(o) \hVVj 



(2.11) 
(2.12) 



Consider the integral in (2.11). In view of (2.10) we can consider locally 



V (z>) = p + -Kz'-z', K = {dtMmjlv 



Therefore we have 



- y(z)| = |p(z') - p\ < ci (/iv^) 2 , z G 5^(0) 



(2.13) 
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which leads to 



< c 2 h VV. 



Suppose |x — hm\ > 2r^h\/V. Since |(x — hm) — z| > |x — hm\/2, we obtain 

(I h ^ W <r ( h ^ H 
K\{(x. — hm) — z) dz < c— 

If |x — hm\ < c lr§h\[V then, for a > r , 



(2.14) 



/ 



x — hm\ n 2 



L 



«a((x — hm) — z) dz 



r h 



< c 



dz 



+ c 



dz 



W-H l(x " hm) ~ Z| ^ 2 \^C)\W-H l(x - ^m) - z 



n-2 ' 



Obviously the first integral in the right-hand side is 0((/i\/^) 2 )- To estimate the second integral 
we use the relation 



|x - z| > c(|x' - z'| + hVv), VzeR n \ # aW ^(x), 
which implies ([3, (4.2.13)]) 

L 



(2.15) 



dz 



x — z 



n-2 



Jb 



dz 



B rohV v{0)\B ahV v(x) (|X' - Z'| + 



n-2 



< c(WP) s 



Thus, for all x 6 M n , 



/ «a((x — /an) — z)dz < c 

Jj 



{hVV) n 



hm\+hVV) n - 2 ' 



Therefore, if Sh denotes the strip S/, = {xe R n : d(x, 50) < rhy/V}, keeping in mind (2.14) 
we have 

L 1 J(,im)l L. . ~ fcm) - z) I" fe) - ' (fi) 



<c\\ f \uhVvr^ Yl 71 — r il,^n-2 



^ch/Hoo^Vp) 



5 h 



5, (|x - y| + /iv 7 ^)"- 2 



^CH/IU^Vp) 2 . (2.16) 



It remains to consider the integral in (2.12). We use the inequality (see [3, p. 80]) 

|k a ((x - ftm) - z) - k a ((x - lira) - y(z)| (2.17) 



9 



< c{hVvf 

and obtain 



(x - hm) - z\ n ~ l + |(x - hm) - y(z)| n " 1 



f |/e A ((x - hm) - z) - « A ((x - /im) - y(z)| |r? (7^=) I dz 



< a (hVv) 2 f -4 i — jdz. 



dz <c (hy/vy 



If |x — hm\ > 2rohy/V then 

If |x — hm\ < 2rohVv, by using (2.15) we obtain that 

f _^ < c (^). 

We deduce that, for all x G W\ 



1 (WP) n 
-ctz < c- 



^o^W l(x-/»m)-z|»-i " (|x-fem| + ^)"- 1 ' 

Therefore 

£ |/(/im)| / |« A ((x - hm) - z)-k a ((x - ftm) - y(z)| |r? f- y(z) < 



<c 1 (hVv) 2 V \f(hm)\[ -4 -— dz 



<c\\f\\ OQ (hVv) 2 f ^— = — 



< C ll/l |oo(/i\/p) 3 | log(max(/i\/p, dist(-x, dQ))\. 
The last inequality and (2.16) complete the proof. 
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3 Implementation and numerical examples 

From Lemma 2.1 and Remark 2.1 we obtain the following one dimensional integral representation 
for the integrals in (2.8) if hm 6 f2 and d(hm, dtt) > rh\/V 



JR" \ h\JV J 



\n/2-i r K n/2 _ 1 {\\x.-y\) (y-hm\ , 

V2M - r— ay 



(2vrP)«/2 J Rn | x - y|n/2-l <™ \ 

h 2 i r (xhVvy 1 / 2 - 1 



= J Rn | rm -y|n/2-l *n/2-l( WD|r m - y|) %M (y) dy 

= i^r e " A2h2 " /4e " irmi2/(i+t)pM(rm ' t)dt 



x — /im 

' ill - 



with 



To apply Lemma 2.1 for nodes hm in the rh\jT> neighborhood of d£l we perform a coordinate 
transformation such that the center hm is the origin in the new coordinates £ = (£i, £ n ) and 
Thm = : £n > Pfem}) where is the distance of the center /im to <9f2 . Since r]2M is radial 
we obtain after the change of variable y — hm = uj£hy/T>, where uj is the rotation matrix, 



X n > 2 ' 1 f tf„/2-i(A|x-y|) /y-/an\ 



(2vr^ ,, ; , m 

{2-KV) n l 2 J( n >^ \U ^m-Cl' 



^" •' e -A 2 h 2 £>t/4 e -|r- m | 2 /(l+t) 



8W 2 7o 

(6rf^ y ■ in'' • / ' •. / •• ~ ' i j i - ■ ■/ i— . •. / v ^ ' in " ' ■ '~™ 



x (erfc (F m (t))P M (a;- 1 r m ,i) + ^ Qm^V^A)) dt, 



where we set 

i J, m(«) = ^(*,(w- 1 r m ) n , P, i 



Then the cubature formula (2.8) reduces to the computation of integrals over the half-line 

,2n r<x> 



^,a/(x) = J] _^ m ) / 00 e-^ a «/*e-l-l a /(i+*)^( rm>t )cft 

d(/im,an)>r-fcv*D ° 

8W 2 ^ ' v Vo 

d(hm,dQ.)<rhVV 

erfc(F m (i))P M (o;- 1 r m ,t) + — ^Q M (^- 1 r m , t, ^)J <ft. 
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For computing at the grid {/ik} we introduce ?7k— m 

= (lo 1 (k — m)) n as the projection of 
k — m onto the normal to the tangential plane, i.e. the £„— axis, and set 



1 J°° e -AV W/4 e -|k|V(D(l+*)) Pm(-^, t) (ft 

i r 

(erfc (F k , m (t))P M (^,t) + ^Qm^ 1 ^,*, ^=)) eft, 



,W = 

4 _ 

b {M) -- / p-A 2 fc 2 77t/4 p -|k-m| 2 /(l?(l+t)) 
k < m " /o 



where 



Hence a cubature for the volume potential at the uniform grid {/ik} is 

Of*) = w4v-r E +^ E ("J 

d(hm,dn)>r hVV d(hm,dQ)<r hy/V 

/im£f2 

We transform the one-dimensional integral representation of d^} m and to integrals over 
M with doubly exponentially decaying integrands by making the substitutions 

t = e ? , £ = a(c7 + e CT ), a = (3(u - e" M ), (3.2) 

with certain positive constants a and b and the computation is based on the classical trapezoidal 
rule. We get after the substitutions 

af > = ly 00 e- A2ft2 ^H/ 4 e-l k l 2 /W+*H))p M (- k =,$(u))$'(n)^ 

fe( M ) = - / e -A 2 fc 2 2?*(«)/4 -|k-m| 2 /(»(l+*(«))) 
k,m g / 

, -F 2 ($(«)) , 

x (erfc(F k , m ($H))^(- ? =-,$(u)) + Qm(u; -i__ *(„), ^L)) ^( u)dUj 

where we set 

3>(u) = exp(a/3(u — exp(— u)) + aexp(/3(-u — exp(— u)))), 
= + e" M )(l + exp(/3(u - exp(-u)))). 

The quadrature with the trapezoidal rule with step size r gives 

iVi , 

E e- A2ft2c *(-)/ 4 e-l k l 2 /^ 1+ *(-)))pM(^^( S r))cI>'( S r), 

s=-AT 

fc( M )«i V e -A 2 ^*(^)/4 e -|k-m| 2 /(C(l+$(sr)))^/ fsrU 
k,m o / y ' V / 

s=-N 

/ k-m p -- F k,m(*( ST )) k-m nu \ 
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We have tested formula (3.1) in the ellipse 



tt = {x€K 2 : 4 + S < 1}» 0<6<a. 
a z b z 

If u(x) = and Vtt(x) = on and outside the ellipse, then the density 

/(x) = (-A + A 2 )u(x), xGfi 

satisfies 

/Ca/(x) 

We have tested the approximation of the potential (1.1) with the density 

/(x) = (-A + A 2 )sin(l-^-|) 2 (3.3) 
which provides the exact value 

/C A /(x) = <j sm ^-^-V?) 



u(x) x G $7 

x $n 



x ^ ft 



and we assumed / itself as the extension outside the ellipse. 

The results in Table 1 show the accuracy of the method. They were obtained with the 
parameters a = 4, (3 = 2 in (3.2) and the quadrature parameters r = 0.01, Nq = —80, Aq = 100. 
We assumed h = 2~ 7 and different a, b and A. We chose r = 6 and V = 3, which gives the 
saturation error less than 10~ 10 . 

In Table 2 and 3 we report on the relative errors and approximation rates for K, \f (0.5, 0) 
and /C\/(0.25, 0.25), respectively, for different a and b, and A = 1. The approximated values are 
computed by the cubature formulas K-^ 1 ^ in (3.1) for M = 1, 2, 3, with P = 4. We have chosen 
a = 4, /3 = 2 in (3.2) and the quadrature parameters r = 0.006, A^o = —160, N\ = 200. We 
obtain rate of convergence 0(h 2M ) although Theorem 2.1. This result comes from the structure 
of the error £ = |/Ca,/ 1 /(x) — /Ca/(x)| which is given by £ = £\ + £2 where £\ = 0(e + h 2M ) and 
£2 = 0(h 2 ). Numerical experiments show that £2 « £\- 

Analogous results are obtained in Table 4 for the approximation of JC\g(0, 0) = 1, with 

g(x u x 2 ) = (-A + A 2 )(l - 4 " 4)V(1 + W 2 )- 1 . (3.4) 

a z b z 

To check the rate of convergence we considered the function 

f(x 1 ,x 2 ) = (-A + A 2 )cos(30 7rxi)cos(307rx 2 ). (3.5) 

The results in Table 5 show that /C^f^/(x) approximates K,\f{~x) with the predicted approxi- 
mation rate 2, for M = 1,2, 3, in agreement with Theorem 2.1. 
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a = b = 


1.5 










X A = 0.2 


A^ = 2 






7 

exact 


approximation 


error 


approximation 


error 


0.00 


0.00 




0.8414709850 


0.258E-09 


0.8414709848 


0.470E-10 


0.25 


0.00 




0.8106234645 


0.267E-09 


0.8106234643 


0.481E-10 


0.50 


0.00 


7104401615 


0.7104401617 


0.300E-09 


0.7104401615 


0.519E-10 


0.75 


0.00 


0.5333026735 


0.5333026737 


0.396E-09 


0.5333026736 


0.687E-10 


1.00 


0.00 


0.3037650630 


0.3037650633 


0.718E-09 


0.3037650631 


0.156E-09 


1.25 


0.00 


0.0932286160 


0.0932286162 


0.248E-08 


0.0932286161 


0.711E-09 


0.25 


0.25 


0.7783135137 


0.7783135139 


0.277E-09 


0.7783135138 


0.492E-10 


0.50 


0.50 


0.5687113034 


0.5687113037 


0.371E-09 


0.5687113035 


0.635E-10 


0.75 


0.75 


0.2474039593 


0.2474039595 


0.895E-09 


0.2474039593 


0.210E-09 


1.00 


1.00 


0.0123453654 


0.0123453656 


0.187E-07 


0.0123453655 


0.581E-08 








a = 1.5, b 


= 1 










A x = 0.2 


A^ = 2 






exact 


approximation 


error 


approximation 


error 


0.00 


0.00 


0.8414709848 


0.8414709867 


0.219E-08 


0.8414709852 


0.519E-09 


0.25 


0.00 


0.8106234643 


0.8106234661 


0.225E-08 


0.8106234647 


0.528E-09 


0.50 


0.00 


0.7104401615 


0.7104401632 


0.248E-08 


0.7104401619 


0.566E-09 


0.75 


0.00 


0.5333026735 


0.5333026752 


0.312E-08 


0.5333026739 


0.679E-09 


1.00 


0.00 


0.3037650630 


0.3037650646 


0.511E-08 


0.3037650634 


0.108E-08 


1.25 


0.00 


0.0932286160 


0.0932286174 


0.154E-07 


0.0932286163 


0.321E-08 


0.25 


0.25 


0.7363058387 


0.7363058405 


0.247E-08 


0.7363058391 


0.580E-09 


0.50 


0.50 


0.3969386023 


0.3969386041 


0.457E-08 


0.3969386028 


0.117E-08 


0.75 


0.75 


0.0351490085 


0.0351490104 


0.541E-07 


0.0351490091 


0.182E-07 








a = 1.5, b 


= 0.5 










A x = 0.2 


A 2 = 2 


XI 


X2 


exact 


approximation 


error 


approximation 


error 


0.00 


0.00 


0.8414709848 


0.8414712254 


0.286E-06 


0.8414710692 


0.100E-06 


0.25 


0.00 


0.8106234643 


0.8106237000 


0.291E-06 


0.8106235459 


0.101E-06 


0.50 


0.00 


0.7104401615 


0.7104403833 


0.312E-06 


0.7104402355 


0.104E-06 


0.75 


0.00 


0.5333026735 


0.5333028743 


0.376E-06 


0.5333027358 


0.117E-06 


1.00 


0.00 


0.3037650630 


0.3037652384 


0.577E-06 


0.3037651117 


0.160E-06 


1.25 


0.00 


0.0932286160 


0.0932287650 


0.160E-05 


0.0932286511 


0.377E-06 


0.25 


0.25 


0.4982722935 


0.4982725309 


0.476E-06 


0.4982723784 


0.170E-06 



Table 1: Exact and approximated values of K,\f{xi,X2) and relative error using IC X h f(x±, X2), 
when h = 2~ 7 and / is given in (3.3) 
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a = 1.5 


b = 1.5 


a = 1.5 


b = 1 


a = 1.5 


b = 0.5 




error 


rate 


error 


rate 


error 


rate 


2 4 


0.439E-01 




0.968E-01 




0.572E+00 




2 5 


0.110E-01 


1.997 


0.243E-01 


1.993 


0.167E+00 


1.781 


2 6 


0.275E-02 


2.000 


0.608E-02 


2.000 


0.419E-01 


1.990 


2 7 


0.688E-03 


2.000 


0.152E-02 


2.000 


0.105E-01 


2.000 


2 8 


0.172E-03 


2.000 


0.380E-03 


2.000 


0.262E-02 


2.000 


2 9 


0.430E-04 


2.000 


0.950E-04 


2.000 


0.655E-03 


2.000 


M = 2 




a = 1.5 


6 = 1.5 


a = 1.5 


b = 1 


a = 1.5 


b = 0.5 


h~ L 


error 


rate 


error 


rate 


error 


rate 


2 4 


0.174E-03 




0.114E-02 




0.186E+00 




2 5 


0.443E-05 


5.294 


0.500E-05 


7.831 


0.288E-02 


6.013 


2 6 


0.183E-06 


4.601 


0.626E-06 


2.996 


0.218E-04 


7.044 


2 7 


0.996E-08 


4.196 


0.534E-07 


3.552 


0.919E-06 


4.570 


2 8 


0.600E-09 


4.053 


0.356E-08 


3.908 


0.922E-07 


3.317 


2 9 


0.371E-10 


4.013 


0.226E-09 


3.979 


0.630E-08 


3.871 


M = 3 




a = 1.5 


b = 1.5 


a = 1.5 


6 = 1 


a = 1.5 


b = 0.5 


ft- 1 


error 


rate 


error 


rate 


error 


rate 


2 4 


0.719E-04 




0.747E-03 




0.469E-01 




2 5 


0.102E-05 


6.138 


0.102E-04 


6.189 


0.177E-02 


4.732 


2 6 


0.155E-07 


6.039 


0.153E-06 


6.066 


0.248E-04 


6.152 


2 7 


0.241E-09 


6.010 


0.236E-08 


6.018 


0.373E-06 


6.057 


2 8 


0.376E-11 


6.001 


0.368E-10 


6.003 


0.577E-08 


6.015 


2 9 


0.936E-13 


5.328 


0.507E-12 


6.180 


0.899E-10 


6.003 



Table 2: Relative errors and approximation rates for /Ca/(0.5,0) = 0.7104401614873481 using 
/CJJV (0-5, 0) in Q = {x e M 2 : § + ff < 1}, A = 1 and / in (3.3) 



M = 1 





a = 1.5 


b = 1.5 


a = 1.5 


6 = 1 


a = 1.5 


6 = 0.5 


h- 1 


error 


rate 


error 


rate 


error 


rate 


2 4 


0.387E-01 




0.955E-01 




0.822E+00 




2 5 


0.967E-02 


2.000 


0.240E-01 


1.991 


0.246E+00 


1.738 


2 6 


0.242E-02 


2.000 


0.601E-02 


1.999 


0.623E-01 


1.983 


2 7 


0.604E-03 


2.000 


0.150E-02 


2.000 


0.156E-01 


1.999 


2 8 


0.151E-03 


2.000 


0.376E-03 


2.000 


0.390E-02 


2.000 


2 9 


0.378E-04 


2.000 


0.939E-04 


2.000 


0.974E-03 


2.000 


M = 2 




a = 1.5 


b = 1.5 


a = 1.5 


6 = 1 


a = 1.5 


6 = 0.5 




error 


rate 


error 


rate 


error 


rate 


2 4 


0.593E-04 




0.139E-02 




0.312E+00 




2 5 


0.225E-05 


4.721 


0.192E-04 


6.174 


0.586E-02 


5.735 


2 6 


0.228E-06 


3.301 


0.244E-06 


6.297 


0.114E-03 


5.682 


2 7 


0.156E-07 


3.870 


0.765E-09 


8.318 


0.350E-05 


5.026 


2 8 


0.997E-09 


3.969 


0.177E-09 


2.113 


0.164E-06 


4.420 


2 9 


0.627E-10 


3.992 


0.145E-10 


3.612 


0.937E-08 


4.126 


M = 3 




a = 1.5 


b = 1.5 


a = 1.5 


6 = 1 


a = 1.5 


6 = 0.5 




error 


rate 


error 


rate 


error 


rate 


2 4 


0.663E-04 




0.762E-03 




0.677E-01 




2 5 


0.947E-06 


6.130 


0.104E-04 


6.191 


0.281E-02 


4.589 


2 6 


0.144E-07 


6.037 


0.156E-06 


6.067 


0.395E-04 


6.154 


2 7 


0.224E-09 


6.009 


0.240E-08 


6.018 


0.593E-06 


6.058 


2 8 


0.348E-11 


6.008 


0.374E-10 


6.004 


0.917E-08 


6.016 


2 9 


0.117E-12 


4.896 


0.558E-12 


6.067 


0.143E-09 


6.003 



Table 3: Relative errors and approximation rates for /Ca/(0.25, 0.25) using K,\ h /(0.25, 0.25) in 

Q = {x G E 2 : § + ff < 1}, A = 1 and / in (3.3). 
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M = 1 





a = 1.5 


b = 1.5 


a = 1.5 


b= 1 


a = 1.5 


b = 0.5 


h- 1 


error 


rate 


error 


rate 


error 


rate 


2 2 


0.415E+00 




0.648E+00 




0.368E+01 




2 3 


0.139E+00 


1.574 


0.216E+00 


1.585 


0.131E+01 


1.489 


2 4 


0.386E-01 


1.854 


0.594E-01 


1.862 


0.366E+00 


1.841 


2 5 


0.992E-02 


1.959 


0.153E-01 


1.961 


0.942E-01 


1.959 


2 6 


0.250E-02 


1.989 


0.384E-02 


1.990 


0.237E-01 


1.989 


2 7 


0.626E-03 


1.997 


0.962E-03 


1.997 


0.594E-02 


1.997 


2 8 


0.157E-03 


1.999 


0.241E-03 


1.999 


0.149E-02 


1.999 


2 9 


0.391E-04 


2.000 


0.602E-04 


2.000 


0.372E-03 


2.000 


M = 2 




a = 1.5 


b = 1.5 


a = 1.5 


b = 1 


a = 1.5 


b = 0.5 


h- L 


error 


rate 


error 


rate 


error 


rate 


2 2 


0.134E+00 




0.210E+00 




0.141E+01 




2 3 


0.173E-01 


2.952 


0.254E-01 


3.050 


0.187E+00 


2.920 


2 4 


0.143E-02 


3.599 


0.207E-02 


3.616 


0.138E-01 


3.756 


2 5 


0.968E-04 


3.882 


0.140E-03 


3.885 


0.914E-03 


3.918 


2 6 


0.618E-05 


3.969 


0.895E-05 


3.970 


0.580E-04 


3.978 


2 7 


0.388E-06 


3.992 


0.562E-06 


3.992 


0.364E-05 


3.995 


2 8 


0.243E-07 


3.998 


0.352E-07 


3.998 


0.228E-06 


3.999 


2 9 


0.152E-08 


3.999 


0.220E-08 


3.999 


0.142E-07 


4.000 



M = 3 





a = 1.5 


b = 1.5 


a = 1.5 


6 = 1 


a = 1.5 


b = 0.5 




error 


rate 


error 


rate 


error 


rate 


2 2 


0.495E-01 




0.703E-01 




0.577E+00 




2 3 


0.284E-02 


4.120 


0.394E-02 


4.157 


0.229E-01 


4.659 


2 4 


0.751E-04 


5.244 


0.106E-03 


5.221 


0.499E-03 


5.517 


2 5 


0.138E-05 


5.768 


0.195E-05 


5.762 


0.903E-05 


5.789 


2 6 


0.225E-07 


5.938 


0.318E-07 


5.937 


0.145E-06 


5.960 


2 7 


0.355E-09 


5.984 


0.502E-09 


5.984 


0.228E-08 


5.990 


2 8 


0.549E-11 


6.014 


0.790E-11 


5.992 


0.357E-10 


5.998 


2 9 


0.150E-12 


5.190 


0.858E-13 


6.524 


0.576E-12 


5.954 



Table 4: Relative errors and approximation rates for lC\g(0, 0) = 1 using JC X h g(0, 0) in Q = 
{x G M 2 : § + ff < 1}, A = 1 and g in (3.4). 
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M = 1 







a = 1.5 6 = 1.5 


a = 1.5 6 = 1 


a = 1.5 6 = 0.5 




xi 


error rate 


error rate 


error rate 


2« 


0.250 
0.750 
1.250 
1.500 


0.171E-05 
0.395E-05 
0.Z68E-U5 
0.569E-08 


0.171E-05 
0.395E-05 

U.Z08E-U0 
0.592E-08 


0.171E-05 
0.395E-05 

U.Z08E-U0 

0.688E-08 




0.250 
0.750 
1.250 
1.500 


0.762E-06 1.163 
0.176E-05 1.163 

n 1 onT7 ncc 11 go 
U.IzUHj-Ud l.loz 

0.140E-07 -1.302 


0.762E-06 1.163 
0.176E-05 1.163 

n 1 omp ncc 11 £;n 

0.144E-07 -1.284 


0.762E-06 1.163 
0.176E-05 1.162 

n 1 omp ncc 1 i go 

0.163E-07 -1.241 


2 iu 


0.250 
0.750 
1.250 
1.500 


0.217E-06 1.811 
0.503E-06 1.811 
0.341E-06 1.811 
0.460E-08 1.609 


0.217E-06 1.811 
0.503E-06 1.811 
0.341E-06 1.811 
0.467E-08 1.627 


0.217E-06 1.811 
0.503E-06 1.811 
0.341E-06 1.811 
0.504E-08 1.690 


2 11 


0.250 
0.750 
1.250 
1.500 


0.561E-07 1.954 
0.130E-06 1.954 
0.881E-07 1.954 
0.122E-08 1.915 


0.561E-07 1.954 
0.130E-06 1.954 
0.881E-07 1.954 
0.123E-08 1.926 


0.561E-07 1.954 
0.130E-06 1.954 
0.881E-07 1.954 
0.128E-08 1.975 


M = 2 






a = 1.5 6=1.5 


a = 1.5 6 = 1 


a = 1.5 6 = 0.5 


ft- 1 


X\ 


error rate 


error rate 


error rate 


2* 


0.250 
0.750 
1.250 
1.500 


0.316E-05 
0.731E-05 
0.497E-05 
0.534E-07 


0.316E-05 
0.731E-05 
0.497E-05 
0.555E-07 


0.316E-05 
0.731E-05 
0.496E-05 
0.644E-07 


2 a 


0.250 
0.750 
1.250 
1.500 


0.897E-06 1.816 
0.208E-05 1.816 
0.141E-05 1.816 
0.194E-07 1.464 


0.897E-06 1.816 
0.208E-05 1.816 
0.141E-05 1.816 
0.196E-07 1.503 


0.897E-06 1.816 
0.208E-05 1.816 
0.141E-05 1.815 
0.210E-07 1.619 


2 iu 


0.250 
0.750 
1.250 
1.500 


0.227E-06 1.986 
0.524E-06 1.986 
0.356E-06 1.986 
0.494E-08 1.970 


0.227E-06 1.986 
0.524E-06 1.986 
0.356E-06 1.986 
0.495E-08 1.984 


0.227E-06 1.986 
0.524E-06 1.986 
0.356E-06 1.986 
0.504E-08 2.056 


2 U 


0.250 
0.750 
1.250 
1.500 


0.567E-07 1.999 
0.131E-06 1.999 
0.890E-07 1.999 
0.124E-08 1.999 


0.567E-07 1.999 
0.131E-06 1.999 
0.890E-07 1.999 
0.124E-08 2.001 


0.567E-07 1.999 
0.131E-06 1.999 
0.890E-07 1.999 
0.124E-08 2.022 


M = 3 






a = 1.5 6=1.5 


a = 1.5 6 = 1 


a = 1.5 6 = 0.5 


ft- 1 


X\ 


error rate 


error rate 


error rate 


2" 


0.250 
0.750 
1.250 
1.500 


0.356E-05 
0.823E-05 
0.559E-05 
0.751E-07 


0.356E-05 
0.823E-05 
0.559E-05 
0.770E-07 


0.355E-05 
0.823E-05 
0.559E-05 
0.862E-07 


2 y 


0.250 
0.750 
1.250 
1.500 


0.906E-06 1.972 
0.210E-05 1.972 
0.142E-05 1.972 
0.198E-07 1.925 


0.906E-06 1.972 
0.210E-05 1.972 
0.142E-05 1.972 
0.198E-07 1.959 


0.906E-06 1.972 
0.210E-05 1.972 
0.142E-05 1.972 
0.203E-07 2.087 


2 1U 


0.250 
0.750 
1.250 
1.500 


0.227E-06 1.999 
0.525E-06 1.999 
0.356E-06 1.999 
0.495E-08 1.999 


0.227E-06 1.999 
0.525E-06 1.999 
0.356E-06 1.999 
0.495E-08 2.002 


0.227E-06 1.999 
0.525E-06 1.999 
0.356E-06 1.999 
0.495E-08 2.035 


2 11 


0.250 
0.750 
1.250 
1.500 


0.567E-07 2.000 
0.131E-06 2.000 
0.891E-07 2.000 
0.124E-08 2.000 


0.567E-07 2.000 
0.131E-06 2.000 
0.891E-07 2.000 
0.124E-08 2.000 


0.567E-07 2.000 
0.131E-06 2.000 
0.891E-07 2.000 
0.124E-08 2.001 



Table 5: Absolute errors and approximation rates for K,\f{x\, 0) using K,\ h f(xi, 0) in Q = {x € 
: f| + § < 1}, with / given in (3.5). 



17 



References 



[1] F. Lanzara, V. Maz'ya and G. Schmidt, On the fast computation of high dimensional 
volume potentials, Math. Comput., 80, 2011, 887-904. 

[2] F. Lanzara, V. Maz'ya and G. Schmidt, Accurate cubature of volume potentials over 
high-dimensional half-spaces, J. Math. Sciences, 173, 2011, 683- 700. 

[3] V.Maz'ya and T. Shaposhnikova, Theory of Sobolev Multipliers, Springer, 2009. 

[4] V. Maz'ya and G. Schmidt, Approximate Approximations, Math. Surveys and Mono- 
graphs, vol. 141, AMS 2007. 

[5] V. Maz'ya, G. Schmidt and W.L. Wendland, On the computation of multi-dimensional 
single layer harmonic potentials via approximate approximations, Calcolo, 40, 2003, 33-53. 

[6] J. Waldvogel, Towards a general error theory of the trapezoidal rule. Approximation and 
Computation 2008. Nis, Serbia, August 25-29, 2008: 
http : //www . math . ethz . ch/~waldvoge/Projects/ integrals .html 



18 



